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Abstract: A quantitative analysis is presented of the dust acoustic wave instability 
driven by the solar and stellar winds. This is a current-less kinetic instability which 
develops in permeating plasmas, i.e.., when one quasi-neutral electron-ion wind plasma 
in its propagation penetrates through another quasi-neutral plasma which contains dust, 
electrons and ions. 

PACS numbers: 96.50. Ci; 96.50.Dj; 52.27.Lw; 52.35.Fp 



During the previous passage of the comet Halley, waves and oscillations in the cometary 
plasma, some of them with very low frequencies in the range of 10 Hz, were clearly 
observed. ^'^Such a cometary plasma itself is a rather complex system which includes 
charged grains and various ion species^, depending on the distance from the coma. Those 
include ions hke H2 , C+, 0+, OH^ , H20^ , etc, with pick-up ions from the solar 
wind being most abundant in the outer regions. The presence of the solar wind plasma, 
which penetrates the cometary cloud, adds to the complexity and makes the whole system 
potentially unstable due to the free energy stored in the streaming solar wind species. 

In our recent work^ a current- less instability of the ion acoustic (lA) perturbations was 
demonstrated in the case of one quasi-neutral electron-ion plasma propagating through 
another static quasi-neutral (target) plasma. The threshold velocity of the propagating 
plasma was shown to be in principle even below the ion acoustic (lA) speed of the static 
plasma, the latter (lA speed) being known to be the threshold speed for the classic case 
of the electron current driven lA wave instability. Hence, the acoustic instability in such 
permeating plasmas may be much more effective than in classic examples with current- 
carrying plasmas, and it may be frequently expected in space and astrophysical plasmas 
where such configurations are rather common. Being current-less, it also removes the 
problem of the self-induced magnetic field associated with the background electron current 
in the classic instability. In the present work this will be investigated in the case of solar 
(stellar) winds propagating through the surrounding interplanetary (interstellar) plasmas 
that also frequently include dust, which may be of cometary or any other origin.^ This 
implies acoustic perturbations which involve the dust dynamics, i.e., the dust acoustic 
mode that has continuously been in the focus of researchers^"^^ ever since its prediction^^ 
and experimental verification.^^. 

One may start from the linearized Vlasov-Boltzmann kinetic equation for the per- 
turbed distribution function 



Here, Ujo — const, and v^^ — nTj/mj, and we assume longitudinal perturbations ~ 
exp{—iuj + ikz) propagating in the direction of the interplanetary magnetic field vector. 
In such a geometry the effects of the magnetic field will not affect the perturbations. 
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where the plasma distribution function for the species j is 
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The perturbed number density may be calculated from 
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For the general species j this yields 



For non-streaming species aj = u/{kvTj), and 



a,- 



rfeexp(-eV2)/(a,-0- 



(2) 



(3) 



(27r)V2 

The integration in is along the Landau contour, and ^ = Vz/v^j- 

For the streaming species the derivation is similar and Eq. ([2]) is obtained, but instead 
of ttj and ^ now we have /3j = {co — kvo) / (kvTj) and C = {vz—vo)/vrpj- The quasi-neutrality 
in the perturbed state 

riwii + n-cii = Uyjei + n^ei + ZdUdi (4) 

will directly yield the dispersion equation. Here, and further in the text, for the cometary 
and wind plasma the indices c and w are used. 

In Eq. ^ for dust the following expansion (to be supported by numbers in the text 
below) will be used 
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This is valid if \ad\ ^ 1 and \Re{a.)d\ ^ \Im{ad)\- For the two electron populations we 
shall use Z{ae) ^ — z(7r/2)^/^ae, ^(/3e) — —iij^/'^Y^'^Pe, that is valid for 

\uj — /cfol 



« 1, |/3e| = 



< 1. 



The same will be used for the cometary ions, \ai\ = \uj\/{kvTci) ^ 1- As for the wind 
ions, the following two interesting cases can be studied. 

The parameter P^ji = {u — kvQ)/{kvTmi) contains two terms, where for the first one we 
expect that uj/{kvTwi) ^ 1, while for the second one in this case we assume fo/^Tu,i ^ 1- 
Hence, the expansion similar to ^ should be used. The dispersion equation in general 
form reads 
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Td \ ^\2J A:.;,,— --^'^ , = i?e[A(fc, a;)] + z/m[A(/c, c)] = 0. (6) 

In application to the solar wind interaction with cometary dusty plasma, one can make 
a few simplifications in Eq. The cometary plasma parameters are used from Ref. 17. 
These imply the following: T^^ = 1.16 ■ 10^ K (= 10 eV), T,i = 2.32 ■ 10^ K (= 2 eV), 



Td = 1.16 ■ 102 K (= 0.01 eV), n^io = 10^ m" 



ndo 



10 m-3, Zd = 800, rud = 1.13 • lO^^o 



kg, and the mean dust grain radius = 3 • 10^^ m. Similar data about the 
dust mass may be found also in Ref. 18 dealing with the latest passage of the comet 
Halley. As for the electrons and ions from the solar wind, we use the following parameters: 
Tew = Tiw = 1.5 ■ 10^ K (~ 13 eV), n^io = fiweo = 5 ■ 10^ m^^. The solar wind speed is 
adopted to be Vweo = "^^uiio = f o = 5 ■ 10^ m/s. Singly charged ions are assumed in both 
systems, and the grains are negatively charged so that 



3.52 ■ 10^ m/s, Vj 



riceo + ZdUdo = ricio- (7) 

The corresponding thermal speeds that will be used 
1.38 ■ 10^ m/s, Vrwe = 1-5 ■ 10^ m/s, v^ce = 1-3 ■ 10^ m/s, Vrd = 0.38 m/s. 

Note that for the interplanetary magnetic field of -Bq = 5 ■ 10^^ T, the corresponding 
gyro-radii for the two populations of protons and electrons are p^i = 74 km, pd = 29 km, 
Pwe = 1-7 km. Pee = 1.5 km. 

With all this it is seen that, first, the contribution of the real term from the ion wind 
part (the terms with the wi-index) is negligible. Second, the ratio of the ion and electron 
imaginary terms is 2;^j(Ti„e/7'w,j)^/2(m^j/me)^/2 exp[— t>Q/(2t>|,^j)] that is completely negli- 
gible. So the total contribution of the wind ions can be omitted. Third, the imaginary 
term from the cometary electrons is negligible as compared to the imaginary term from the 
cometary ions; the ratio of the two is proportional to {nceo/ncio){Tei/TceY^'^{me/'meiy^'^ ^ 
1. 

As a result, the real part of ([6]) yields the frequency of the dust acoustic mode 
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The growth rate 7 ~ —ImA{k,Ur)/[d{ReA)/duj]ujc^ujr is 



7 



where 



T 



1 + 



^weO 



TT 



1/2 



3/2 



1/2 



me, 
rrip 



1/2 



(9) 



4 



Hence, the instability sets in if 



vo > -^(1 + a), 
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For the parameters given above the threshold velocity Vth for the instability is in fact very 
low. Taking Zct = 800, from (fTOj) we have Vth = 5.3 km/s only. Therefore, within the 
present model the wind-driven dust acoustic oscillations are always growing. In Fig. 1 
the growth rate is given in Hz, in terms of the wave number and the dust charge number 
density. 

The wave frequency contour plot has the same shape, with the maximum frequency 
around 2 Hz. In the case of the fast solar wind vq = 8 ■ 10^ m/s and for Zd = 800, the 
ratio ■y/ur is increased by about 60 percent. 

In such a multi-component system, the Debye length is determined by the coolest 
species (the dust in the present case), and it turns out to be of the meter size. Hence, 
the shortest wavelengths of interest here should be above that scale and consequently the 
largest wave frequency that should be attributed to the DA mode is expected to be of the 
order of a few Hz. In view of large gyro-radii given previously in the text, the magnetic 
field effect is clearly completely negligible even for the lightest species. 

In the three-component cometary plasma the DA phase speed is 
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For the same parameters as above we have c^^ = 3.9 m/s. In the presence of the wind, 
the actual DA phase speed Vph = Ur/k = 3.8 m/s. Observe that the threshold velocity is 

Vth = ^-^Vph- 

In the case 

\u - kvo\/{kvT^,) < 1 
the expansion for the wind ions is different and it yields 
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A procedure similar as above yields the frequency 
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yields the instability condition which in this case reads 
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For the solar wind data the parameter 6 ~ 43, so the critical velocity ( !T5|) is reduced. 

To conclude, the demonstrated current-less instability is expected to be very effective 
in the environment of cometary tails, and in interstellar dusty plasma penetrated by 
the stelar winds in general. The threshold velocity is very low, practically equal to the 
phase speed of the DA wave propagating in such a two-plasma system [this particularly in 
the limit f|TT]) ]. and the restoring pressure appears mainly due to the contribution of the 
cometary electrons. The model presented here is based on the assumption of a fixed charge 
on dust grains. However, in realistic situations inelastic collisions (i.e., those related to 
charging of dust grain and absorption of plasma particles) may play an important role and 
yield considerably different results. These effects have been studied in detail in Ref. 13 
and in several references cited therein. In the range of large wave-lengths the damping 
rate due to dust charging can become several orders of magnitude larger than the Landau 
damping on ions, while in the short wave-length domain the charging effects are almost 
absent. This short wave-length range is in fact of interest in the present study because 
it provides a much more strongly growing mode, as can be seen from Fig. 1, and in this 
regime the charging phenomena should not considerably alter the mode behavior. 



6 



References 

[I] R. Grard, A. Pedersen, J. G. Trotignon, C. Beghin, M. Mogilevsky, Y. Mikhailov, O. 
Molchanov, and V. Formisano, Nature 321, 290 (1986). 

[2] S. Klimov, S. Savin, Ya. Aleksevich, G. Avanesova, V. Balcbanov, M. Balikhin, A. 
Galeev, B. Gribov, M. Nozdrachev, V. Smirnov, A. Sokolov, O. Vaisberg, P. Oberc, Z. 
Krawczyk, S. Grzedzielski, J. Juchnicwicz, K. Nowak, D. Orlowski, B. Parfianovich, D. 
Wozniak, Z. Zbyszynski, Ya. Voita, and P. Triska, Nature 321, 292 (1986). 

[3] H. Balsiger, K. Altwcgg, F. Bhlcr, J. Geiss, A. G. Ghielmctti, B. E. Goldstein, R. 
Goldstein, W. T. Huntress, W. H. Ip, A. J. Lazarus, A. Meier, M. Neugebauer, U. 
Rettenmund, H. Rosenbauer, R. Schwenn, R. D. Sharp, E. G. Shelley, E. Ungstrup, 
and D. T. Young, Nature 321, 330 (1986). 

[4] J. Vranjes, S. Poedts, and Z. Ehsan, Phys. Plasmas 16, 074501 (2009). 

[5] H. Saleem, Phys. Plasmas 13, 012903 (2006). 

[6] C. B. Dwivedi and B. P. Pandey, Phys. Plasmas 2, 4134 (1995). 

[7] J. Vranjes, B. P. Pandey, and S. Poedts, Phys. Rev. E 64, 06640416 (2001). 

[8] V. E. Fortov, A. D. Usachev, A. V. Zobnin, V. I. Molotkov, and O. F. Petrov, Phys. 
Plasmas 10, 1199 (2003). 

[9] S. Mahmood and H. Saleem, Phys. Lett. A 338, 345 (2005). 

[10] J. Vranjes and S. Poedts, Phys. Plasmas 13, 052103 (2006). 

[II] T. K. Baluku and M. A. Hellberg, Phys. Plasmas 15, 123705 (2008). 
[12] J. Vranjes and S. Poedts, Phys. Rev. E 82, 026411 (2010). 

[13] P. Tohas S. Ratynskaia, and U. de Angehs, Phys. Plasmas 17, 103707 (2010). 
[14] M. Tribeche and A. Merriche, Phys. Plasmas 18, 034502 (2011). 
[15] N. N. Rao, P. K. Shukla, and M. Y. Yu, Planet. Space Sci. 38, 543 (1990). 
[16] A. Barkan, N. D'Angelo, and R. L. Merhno, Planet. Space Sci. 44, 239 (1996). 
[17] N. D'Angelo, Planet. Space Sci. 46, 1671 (1998). 

7 



[18] J. A. M. McDonnell, W. M. Alexander, W. M. Burton, E. Bussoletti, D. H. Clark, 
R. J. L. Grard, E. Grn, M. S. Manner, D. W. Hughes, E. Igenbergs, H. Kuczera, B. A. 
Lindblad, J. C. Mandeville, A. Minafra, G. H. Schwehm, Z. Sekanina, M. K. Wallis, 
J. C. Zarnecki, S. C. Chakaveh, G. C. Evans, S. T. Evans, J. G. Firth, A. N. Littler, 
L. Massonne, R. E. Olearczyk, G. S. Pankiewicz, T. J. Stevenson, and R. F. Turner, 
Nature 321, 338 (1986). 

[19] V. N. Tsytovich and U. de Angelis, Phys. Plasmas 6, 1093 (1999); ibid. 
7, 554 (2000); ibid. 8, 1141 (2001). 



8 




0.1 0.2 0.3 0.4 

k [1/m] 

Figure 1: The growth rate (in Hz) of the DA mode in a cometary dusty plasma, driven 
by the solar wind, in terms of the wave-number and the grain charge number. 
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